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ABSTRACT 


An  analytical  model  of  a  jet  injected  normally  from  a 
flat  plate  into  a  uniform  crossing  flow  was  modified  to 
provide  an  improved  method  of  predicting  the  interference 
effects  arising  from  the  complex  flow  fields  induced  by  ship 
bow-thrusters.  This  model  is  an  extension  of  previous  work 
based  upon  a  description  of  the  jet  as  a  series  of  distrib¬ 
uted  vortices.  The  analysis  takes  into  account  the  position 
of  the  effective  source  of  the  jet  and  the  blockage  due  to 
the  presence  of  the  jet  in  the  crossflow.  Improvement  in 
the  predicted  pressure  distribution  on  the  plate  was  obtain¬ 
ed  by  adopting  a  new  formula  for  determining  the  effective 
jet  source  which  more  closely  approximates  the  experimental 
results.  Further  improvement  resulted  from  the  incorpora¬ 
tion  of  a  symmetric  foil  shape  to  simulate  the  presence  of 
the  jet  and  its  highly  turbulent  wake.  The  accuracy  of  the 
model  was  evaluated  for  representative  jet-to-crossflow 
velocity  ratios.  Good  agreement  with  experimental  results 
was  achieved  for  large  portions  of  the  interaction  field 
with  particular  improvement  in  the  near-field  and  down¬ 
stream  regions. 
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I .  INTRODUCTION 


Since  man  first  set  sail,  ships  and  boats  have  been 
steered  by  stern-mounted  rudders  or  rudder-like  devices. 

Die  bow  has  been  left  relatively  uncontrolled.  As  the  size 
of  ships  increased,  they  became  more  difficult  to  maneuver 
in  close  quarters,  as  when  entering  port.  Tug  boats  were 
often  required  to  provide  lateral  thrust  at  both  ends  of  the 
ship.  This  was  a  satisfactory  solution  to  the  problem  for 
many  years.  It  did  not  provide  the  degree  of  independence 
desired  by  ship  operators,  however,  and  other  methods  were 
considered.  Although  experimental  lateral  thrust  units  were 
used  in  the  mid-nineteenth  century,  the  first  practical  unit 
was  not  installed  until  1955.  This  unit,  a  through- tunnel 
type  bow  thruster  installed  in  the  "Princess  of  Vancouver" 
proved  highly  successful  [1] .  Since  then,  bow  thrusters 
have  steadily  gained  in  popularity.  The  advantages  of  these 
units  include  providing  a  greater  degree  of  maneuverability, 
increased  safety,  increased  ability  to  operate  in  rough 
weather  at  oil  rigs  and  research  sites,  and  reduced  tug 
charges. 

In  the  past  twenty-five  years,  ports  have  become  increas 
ingly  crowded,  labor  and  fuel  costs  have  greatly  increased 
tug  boat  charges,  and  ships  have  been  built  that  are  much 
larger  and  less  maneuverable.  Thus  lateral  thrust  units 


have  not  only  become  economical,  but  they  are  standard 
equipment  for  many  new-construction  ships. 

As  more  bow  thrusters  came  into  use,  it  became  apparent 
that  although  they  work  well  at  zero  or  very  low  forward 
speed,  their  effectiveness  drops  considerably  as  the  ship's 
speed  is  increased  [2].  In  view  of  the  large  number  of 
these  units  currently  installed,  there  is  suprisingly  little 
information  available  concerning  thruster/hull  interactions. 

Soon  after  the  first  unit  was  installed  in  1955,  it  was 
learned  that  the  operation  of  a  tunnel-type  thrust  unit 
causes  a  low  pressure  area  to  form  on  the  ship's  hull  [2] 
(see  Fig.  1).  To  find  the  net  resultant  thrust,  this  suc¬ 
tion  force  must  be  integrated  over  the  effected  area  and 
subtracted  from  the  thrust  of  the  unit.  When  a  thrust  unit 
is  operated,  water  surrounding  the  jet  is  entrained,  as  in 
an  eductor,  causing  a  low  pressure  area  to  be  formed  on  the 
ship's  hull  in  a  small  area  about  the  jet.  In  addition,  if 
a  ship  is  moving  forward,  a  plume  is  developed  aft  of  the 
thrust  unit  and  a  system  of  vortices  is  formed.  These  vor¬ 
tices  induce  a  velocity  in  the  water  at  the  ship's  hull, 
causing  a  pressure  drop.  Additionally,  the  jet  itself  acts 
as  an  obstruction  to  the  flow  of  water  past  the  ship,  caus¬ 
ing  the  flow  to  change  speed  as  it  passes,  thus  causing 
another  pressure  variation.  These  induced  velocity  and 
blockage  effects  are  so  severe  that  a  thrust  unit's 


effectiveness  can  drop  as  much  as  fifty  percent  with  a  speed 

increase  of  only  three  knots  [1] . 

The  effects  of  bow  thruster  jet/hull  interactions  upon 

a  ship  motion  can  be  better  understood  by  considering  the 

equations  of  motion  for  a  ship  moving  in  a  horizontal  plane. 

The  general  forms  of  the  equations  are  [3] : 

m(u-vtp)  =  X  Surge  equation 

m(v+uip)  *  Y  Sway  equation 

I  ji  *  N  Yaw  equation 

z 

The  linearized  forms  of  these  equations  are  sufficient 
for  this  analysis  and,  making  the  ship's  center  of  gravity 
the  coordinate  system  reference  point,  they  become 
respectively: 

-Xu(u-uQ)  +  (m-x^Ju  »  X  (1) 

-Yvv  +  (m-Y^)v  -  (Yr-muo)r  -  Y^r  «  Y  (2) 

-Nyv  -  N^v  -  Nrr  +  (Iz-Ni.)r  =  N  (3) 

where : 

m  =  mass  of  the  ship 

u  *  velocity  in  the  x-direction  (initial  value  uQ) 
v  *  velocity  in  the  y-direction 
*  yaw  angle 

r  **  angular  velocity  =  ^ 

I  *  mass  moment  of  inertia 
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X  *  excitation  forces  in  the  x-direction 

Y  »  excitation  forces  in  the  y-direction 

N  *  excitation  torque 

Dots  indicate  differentiation  with  respect  to  time. 

The  subscripts  u,  u,  v,  v,  r,  r,  denote  differenti¬ 
ation  with  respect  to  the  subscript  variable. 

The  left  sides  of  these  equations  contain  the  velocity 
and  acceleration-dependent  terms  such  as  added  mass,  added 
inertia,  and  hydrodynamic  drag.  The  right  sides  implicitly 
represent  the  total  forces  in  the  x  and  y  directions  and  the 
total  moment  about  a  vertical  axis  through  the  center  of 
gravity.  Although  these  terms  are  used  to  represent  various 
excitation  forces  acting  on  a  ship  due  to  rudders,  stabi¬ 
lizers,  propellors,  etc.,  this  analysis  will  deal  solely 
with  the  forces  and  moments  caused  by  a  through-tunnel 
bow-mounted  thrust  unit. 

In  the  y-direction,  the  force  due  to  a  bow  thruster  may 
be  written  as: 

Y  *  Y.  +  y„ 

t  s 

where : 

Y  a  -m.u. ,  the  ideal  jet  thrust 

^  J  J 

Y  *  the  induced  "suction" force  on  the  hull 
s 

ihj  *  the  mass  flow  rate  of  the  jet 

Uj  »  the  jet  velocity  (positive  to  starboard) . 
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Inspection  of  Fig.  2  reveals  that  -Y  /Yt  nearly  linear 
in  the  (u^/Ujl  ratio  range  of  .1  to  .3.  Using  this  velocity 
ratio  range  as  an  example,  an  approximation  of  the  right 
side  of  equations  (2)  and  (3)  cam  be  made.  This  portion  of 
the  curves  is  of  interest  because  the  largest  decrease  in 
thruster  performance  tadces  place  here.  Since  the  slope  and 
intercept  of  this  line  are  dependent  on  hull  shape,  am 
average  performance  curve  has  been  plotted  in  Fig.  3  to 
represent  a  "typical"  hull.  A  straight  line  having  a  slope 
of  5.3  amd  am  intercept  of  -0.65  is  used  to  approximate  the 
curve  in  the  I u  /u . I  ratio  (hereafter  R)  of  . 1  to  . 3 . 

Thus: 

-Y/Y.  *  5.3R  -  .65 
s  t 

Then  since 

Yt  -  -mjuj 

Ys  -  *^<5.311  -.65)  (4) 

and 

Y  *  -m_.UjU.65  -5.3R)  .1  a  R  a  .3  (5) 

By  substituting  this  equation  into  equation  (2) ,  the 
motion  in  the  y-direction  is  explicity  coupled  to  the  motion 
in  the  x-direction  through  the  dependency  upon  R. 

The  total  turning  moment  N  is  the  moment  produced  by 
the  force  Y  acting  at  a  distance  x  from  the  center  of  gravity 


This  characteristic  moment  arm  varies  non- linearly  as  the 
ship's  speed  changes  and  is  given  as: 


Vt  ♦  Vs 


where: 

■  the  distance  from  the  center  of  gravity 
to  the  thrust  unit 

xg  *  the  distance  from  the  center  of  gravity 
to  the  center  of  action  of  the  hull 
suction  force. 


Chislett  and  Bjorheden  [2]  found  that  xg  varies  linearly 
with  ship's  speed.  The  curve  in  Fig.  4,  which  is  adapted 
from  their  data,  has  a  slope  of  -.7  and  an  intercept  of  1.0. 
Therefore : 

xg  «  xfc(l  -  . 7R)  R  <  1.2  (6) 

When  equations  (4)  and  (5)  are  substituted  for  Y  ,  Y.  and  Y, 

9  L 

the  location  of  the  center  of  action  of  the  combined  force 
system  is  given  by: 


„  ,  (3.71RZ  -  5.76R  +  1.65),. 

x  x. 

G  (1.65  -  5.3R'  z 


.1  <  R  <  .3 


(7) 


and 


Thus  the  yaw  equation  is  coupled  to  the  sway  and  surge 
equations.  Note  that  equation  (7)  is  only  valid  for  velocity 
ratios  between  .1  and  .3  for  this  particular  fictitious  hull. 

Although  this  analysis  has  been  greatly  simplified,  it 
demonstrates  the  complication  of  the  relationships  between 
turning  moments  and  ship's  speed  when  a  lateral  thrust  unit 
is  used.  It  can  be  seen  from  Fig.  2  that  model  tests  of 
each  ship  type  are  necessary  in  order  to  accurately  predict 
the  performance  characteristics  of  a  lateral  thrust  unit. 

With  a  thorough  understanding  of  the  fluid  mechanics  of  bow 
thruster  systems,  thrust  units  could  be  designed  to  reduce 
the  hull  suction  effect  and  thus  make  them  more  usable  under 
a  wider  range  of  conditions.  An  accurate  analytical  method 
for  predicting  jet /hull  interactions  would  be  useful  in 
supplying  this  information  without  lengthy  and  costly  model 
tests.  Unfortunately,  an  accurate  model  that  can  be  super¬ 
posed  with  other  effects  is  not  available.  It  is  the  purpose 
of  this  thesis  to  propose  an  analytical  method  for  predicting 
the  jet/hull  interactions  of  a  through- tunnel  type  lateral 
thrust  unit. 
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II.  BACKGROUND 


As  demonstrated  in  the  proceeding  section,  the  ship's 
hull/thruster  interactions,  when  a  ship  has  way  on,  cam 
seriously  degrade  a  bow-thruster's  effectiveness.  In  order 
for  naval  architects  to  design  a  bow  thruster  installation 
so  that  it  can  be  used  effectively  when  the  ship  is  moving, 
a  complete  knowledge  of  these  effects  is  required.  Unfortu¬ 
nately,  only  a  few  waterborne  tests  have  been  performed. 

These  experiments  (Fig.  2) ,  while  clearly  demonstrating  the 
overall  effects  on  turning  ability,  do  not  provide  informa¬ 
tion  on  the  actual  pressure  distribution  on  the  hull.  For 

an  incompressible  flow,  the  jet  path  and  the  induced  flow 

2  2 

are  dependent  mainly  on  the  momentum  flux  ratio  p^U.,  /p^ 

and  are  relatively  independent  of  the  Reynolds  number  [4]. 
Thus  an  air  jet  injected  into  a  crossing  air  flow  will  induce 
a  reaction  similar  to  that  of  the  water  jet  of  a  thrust  unit 
at  the  same  velocity  ratio.  Uiere  have  been  a  number  of 
studies  of  air  jets  emitted  from  a  flat  plate  into  a  crossing 
air  flow  in  connection  with  the  design  of  VSTOL  (Vertical  and 
Short  Take  Off  and  landing)  aircraft.  Notable  studies  have 
been  made  by  Bradbury  and  Woods  [4]  and  Feam  and  Weston  [5]  . 

While  the  results  of  these  experimental  studies  are  use¬ 
ful  to  the  naval  architect  when  making  an  initial  design 
intended  to  improve  thruster  efficiency,  they  provide  no 
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method  for  evaluating  the  resultant  configuration.  Lengthy 
and  expensive  model  or  shipboard  tests  are  required  for  this 
purpose.  An  accurate  and  versatile  analytical  model  of  the 
flow  interaction  effects  would  provide  a  quick  and  inexpen** 
sive  means  of  evaluating  new  designs,  thereby  reducing  the 
required  number  of  model  tests.  Some  empirical  and  semi- 
empirical  models  have  been  developed.  While  they  can  predict 
the  pressure  distributions  on  a  plate,  they  cannot  be  easily 
superposed  with  other  effects  such  as  a  secondary  jet. 

Because  a  fluid  jet  injected  into  a  crossing  stream  is  a 
turbulent,  three  dimensional,  highly  non-linear  problem, 
empiricism  to  some  degree  has  been  used  in  all  models  to 
date.  In  this  study  an  effort  has  been  made  to  eliminate  or 
at  least  reduce  the  empirical  inputs  to  existing  analytical 
models  while  improving  their  accuracy. 

A  review  of  available  models  revealed  that  Wooler's  formu¬ 
lation  [61  minimized  empiricism,  requiring  only  an  experi¬ 
mentally  determined  jet  trajectory.  Hence  his  model,  as 
modified  by  Waterman  [7] ,  was  selected  as  a  basis  for  further 
investigation . 

In  order  to  visualize  the  interaction  field  and  determine 
the  accuracy  of  an  improved  model,  once  formulated,  accurate 
experimental  results  were  required.  A  review  of  available 
studies  revealed  that  the  pressure  contours*  for  a  given 

* 

In  this  thesis,  values  of  local  pressures,  P,  are  expres¬ 
sed  in  the  form  of  pressure  coefficients,  c  ,  where: 

P-P  p 


velocity  ratio  vary  considerably  from  one  report  to  the  next, 
as  illustrated  in  Fig.  5  and  6.  This  difference  is  especially 
noticeable  in  the  downstream  region  for  values  of  |n^/Uaj 
(hereafter  VR)  of  about  4 ,  and  is  present  to  some  degree  at 
all  velocity  ratios.  Figure  5  seems  to  Imply  that  the  induced 
flow  may  not  be  solely  a  function  of  velocity  ratio  since  the 
different  studies  were  not  made  at  the  same  freestream 
velocities.  However,  considering  studies  made  with  the  cross- 
flow  velocity  nearly  equal,  as  in  Fig.  6,  it  is  seen  that  the 
discrepancies  are  more  probably  due  to  variations  in  the 
equipment  and  test  procedures  used.  This  serves  to  point  out 
the  uncertainty  of  the  available  experimental  results.  The 
results  of  Feam  and  Weston  [5]  agree  well  with  those  of 
other  researchers  over  a  broad  range  of  velocity  ratios. 
Additionally  their  report  was  the  most  complete  of  those 
reviewed.  It  contains  not  only  the  usual  plots  of  pressure 
contours,  but  also  raw  data,  time  variation  of  pressure 
information,  and  a  comparison  with  other  reports.  For  these 
reasons,  the  results  reported  in  [5]  were  selected  as  the 
basis  for  determining  the  accuracy  of  the  analytical  model. 

Before  attempting  to  develop  an  improved  model,  a  thor¬ 
ough  review  of  Wooler's  and  Waterman's  work  was  made  in 
order  to  understand  their  approaches  and  verify  their  results. 
A  rederivation  of  these  basic  formulations  will  not  be  pre¬ 
sented  here  as  it  can  be  found  in  references  [6]  and  [7] .  It 
was  noted  that  Waterman's  adoption  of  a  conforming  vortex 
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model  was  a  considerable  improvement  to  Wooler's  formulation. 
In  the  conforming  vortex  model  the  vortices  are  constrained 
to  lie  within  the  jet  rather  them  being  tangential  to  the  jet 
axis  as  in  wooler's  model.  By  making  this  alteration,  the 
horseshoe  vortices  are  made  to  conform  to  the  curvature  of 
the  jet  described  by  the  empirical  trajectory  equation  (see 
Pig.  7)  : 


|  -  B [cosh (Z/Bd)  -  1]  (8) 

Where: 

x  *  the  distance  down-stream  from  the  jet  center. 

d  *  the  jet  diameter. 

z  ■  height  of  the  jet  axis  above  the  plate. 

2 

B  *  trajectory  parameter  »  .  19(VR)  . 

The  interference  velocity  at  an  arbitrary  point  on  the 
plate  is  determined  by  applying  the  Biot-Savart  Law.  The 
trailing  vortices  are  approximated  as  straight  line  segments. 
As  these  segments  become  small,  the  linear  approximation  can 
be  made  to  conform  to  the  curvature  of  the  jet.  The  total 
interference  component  at  a  particular  point  due  to  the 
entire  jet  results  from  the  summation  of  all  the  elemental 
contributions . 

A  comparison  of  the  results  of  both  formulations  reveals 
that  Wooler's  work  more  closely  agrees  with  the  experimental 
findings  of  Bradbury  and  Woods  [4]  them  does  Waterman's, 
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even  though  the  trailing  vortices  in  Wooler's  model  do  not 
follow  the  jet  trajectory  as  in  the  actual  case.  This  result 
is,  in  part,  due  to  the  fact  that  in  wooler's  formulation  the 
perturbation  velocity  at  a  particular  point  on  the  plate  is 
caused  by  vortices  acting  at  different  distances  from  the 
plate.  In  Waterman's  model  the  interference  velocities 
result  from  the  trailing  horseshoe  vortex  contributions  act¬ 
ing  at  a  single  radius.  Thus  the  interference  velocities 
predicted  by  Waterman's  conforming  vortex  model  are  generally 
higher  and  hence  produce  lower  coefficients  of  pressure  than 
those  of  Wooler.  In  addition  to  modifying  wooler's  model, 
and  noting  studies  by  Keffer  and  Baines  [8]  and  Schmidt  [9] , 
Waterman  introduced  the  concept  of  an  effective  origin  of  the 
vortex  sheet  at  a  distance,  Z  ,  above  the  plate  given  by 

w 

the  furmula: 

„  _  VR  VW 

ze  — rr~  (9) 

Applying  this  expression  to  the  conforming  vortex  model. 
Waterman  produced  the  same,  if  not  better  far-field  agreement 
(Pig.  8).  However,  the  model  lacked  accuracy  in  the  near¬ 
field. 

Waterman  further  noted  that  the  discrepancies  observed  in 
the  upstream  region  can  be  at  least  partially  attributed  to 
the  blockage  of  the  freest ream  flow  by  the  jet.  It  was  noted 
that  the  downstream  area  is  dominated  by  separation,  vortex 
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III.  IMPROVEMENT  IN  THE  EXISTING  MODEL 


Evaluation  of  Waterman's  "fully-corrected  conforming  vor¬ 
tex  model"  (Figs.  9  and  10)  revealed  that  the  correlation  of 
the  analytical  model  with  the  experimental  data  in  the  area 
just  forward  of  the  jet  is  adequate  at  each  of  the  velocity 
ratios  reported.  However,  his  model  is  accurate  in  the  far- 
field  only  for  the  velocity  ratio  of  8.0  and  in  the  near¬ 
field  to  the  sides  of  the  jet  for  a  velocity  ratio  of  4.0. 

The  correlation  is  poor  in  the  area  aft  of  the  jet  for  all 
velocity  ratios.  It  was  apparent  that  the  near  and  mid- 
field  regions,  especially  aft  of  the  jet,  required  the 
greatest  attention. 

In  order  to  gain  a  better  understanding  of  the  effects  of 
the  empirical  inputs  to  Waterman's  model,  it  was  run  a  num¬ 
ber  of  times  using  various  values  of  the  constants  in  the 
trajectory,  effective  source,  and  effective  diameter 
equations.  It  was  learned  that  small  changes  in  the  trajec¬ 
tory  parameter,  B  in  equation  (8),  cause  little  change  in 
the  overall  pressure  distribution.  The  correlation  with 
experimental  data  was  improved  at  some  velocity  ratios  while 

it  diverged  at  others  for  a  particular  value  of  B.  A  value 
2 

of  B  ■  . 19 (VR)  ,  as  before,  resulted  in  the  best  overall 
correlation  of  the  pressure  distribution  on  the  plate. 
Changing  the  effective  diameter  of  the  jet  causes 


22 


i 

f 

i 


major  changes  in  the  near-field  pressure  distribution,  but 
has  little  effect  in  the  far-field.  This  was  expected  con¬ 
sidering  the  short  range  of  the  disturbance  caused  by  a 
doublet  in  a  uniform  flow  using  inviscid  relationships.  The 
model  proved  to  be  particularly  sensitive  in  both  the  near 
and  far-field  regions  to  changes  in  the  position  of  the 
effective  source.  Since  the  quality  of  the  model  depends 
so  greatly  on  the  effective  source  correction,  the  equation 
for  it  was  investigated  more  closely. 

Equation  (9) ,  the  empirical  formula  for  the  effective 
source  height,  was  formulated  by  Schmidt  [9]  who  refers  to 
Jordinson  [10]  and  Keffer  and  Baines  [8]  for  the  experimental 
data  used  as  a  basis  for  his  work.  The  Jordinson  paper  could 
not  be  obtained  for  review.  The  data  presented  by  Keffer 
and  Baines,  however,  does  not  fit  the  equation  presented  by 
Schmidt.  The  equation,  Z  /d  *  . 8 "V  VR  fits  the  Keffer  and 
Baines  data  much  better  as  seen  in  Table  1. 


TABLE  I 

COMPARISON  OF  VALUES  OF  Zfi/d 


VR 

Schmidt 

Experimental 
Keffer  & 
Baines 

.  8  VW 

4 

.53 

1.59 

1.60 

6 

.97 

1.95 

1.93 

8 

1.51 

2.30 

2.26 
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However,  Z  /d  can  vary  from  .  54  v  VR  to  .8VWand  still  be 

in  the  range  of  the  heights  given  in  Table  I.  Using  an 

ef f ective  source  location  of  Z  /d  »  .8  *  VR  in  the  Waterman 

e 

model  significantly  improves  the  correlation.  But,  if 
Z  /d  ■  .6  v  is  used,  the  correlation  is  much  better  even 
without  the  effective  diameter  correction.  Thus,  by  estab¬ 
lishing  the  effective  source  at  . 6  V VR  diameters  above  the 
plate,  the  full  jet  diameter  can  be  used  in  the  model  and 
the  diameter  correction  is  eliminated.  See  Fig.  11  for  an 
example  of  the  improved  correlation  resulting  from  this 
change . 

Concentrating  on  the  area  aft  of  the  jet,  it  can  be  seen 

by  comparing  Fig.  8  with  Fig.  10  that  this  area  is  primarily 

affected  by  the  blockage  near  the  plate  due  to  the  jet.  In 

making  multiple  runs  of  Waterman’s  model  and  varying  the 

parameters  of  the  vortices  within  the  plume,  it  was  confirmed 

that  they  have  little  effect  on  the  near-field  pressures. 

The  far-field  to  either  side  of  the  jet  is  dominated  by  the 

vortices  within  the  plume.  Therefore,  to  improve  the  model's 

accuracy  in  the  near-field,  the  analysis  centered  on  the  jet 

blockage  effect.  Wu  and  Wright  [11]  proposed  that  a 

fictitious  afterbody  exists  downstream  of  the  jet.  This 

afterbody  encompasses  the  highly  turbulent,  three-dimensional 

flow  caused  by  vortex  shedding  and  entrainment  of  fluid 

within  and  downstream  of  the  jet.  Outside  of  this  area,  the 

fluid  streamlines  are  fairly  smooth  and  close  to  those  of  an 

inviscid  flow.  With  this  concept  in  mind,  consideration  was 
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given  to  using  an  elongated  shape  such  as  an  ellipse  or  a 
symmetrical  foil  to  approximate  the  blockage  effect. 

In  order  to  establish  a  probable  shape  for  the  afterbody, 
the  coefficients  of  pressure  along  the  corssflow  axis  from 
the  forward  edge  of  the  jet  to  seven  diameters  aft  of  the 
jet  were  plotted  using  the  data  of  [5]  (Fig.  12) .  These 
pressure  contours  resemble  those  in  reference  [12]  and  [13] 
for  a  symmetrical  foil  in  a  uniform  flow.  Thus  a  symmetrical 
foil  was  selected  to  approximate  the  jet  and  the  afterbody 
near  the  plate. 

To  determine  a  length  for  the  foil,  it  was  assumed  that 
the  foil  should  end  where  the  coefficient  of  pressure  reaches 
an  undisturbed  (zero)  value.  It  was  not  expected  that  the 
coefficient  of  pressure  should  necessarily  rise  above  zero 
aft  of  the  jet  due  to  severe  dissipation  in  the  near-wake 
region  and  the  interference  effects  of  the  vortices  within 
the  plume.  Indeed  the  coefficient  of  pressure  leveled-off 
at  approximately  -0.1  for  all  available  data.  Therefore, 

-0.1  was  selected  as  a  first  approximation  of  the  value  of 
the  coefficient  of  pressure  used  to  determine  the  length  of 
the  foil.  The  corresponding  distances  from  the  jet  center 
to  the  trailing  edge  were  determined  to  be  6,2,  3.0,  and  2 . 5 
jet  diameters  for  the  velocity  ratios  of  3.9,  8.Q,  and  10.0, 
respectively. 

For  simplicity,  a  standard  Joukowsky  foil  was  selected 
since  it  can  be  easily  approximated  through  mapping  methods. 
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Referring  to  Fig.  13,  the  Joukowsky  transformation. 


Z  »  c 


is  used  to  transform  the  foil  in  the  z-plane  to  a  cylinder 
in  the  c-plane  so  that  the  foil  can  be  treated  as  a  doublet 
in  a  uniform  stream.  The  perturbation  velocities  due  to  the 
presence  of  the  foil  are  then: 


u 

U 

OB 


(c'  +  a) 
(C  +  i)~ 


1 


and 


2 


(c'  +  a) 
(c  +  D 


Where: 

u  =  the  surface  perturbation  velocity  in  the  direction 
of  the  cross- flow  due  to  the  presence  of  the  jet 
and  afterbody  as  modeled  by  the  foil. 

v  »  the  surface  perturbation  velocity  perpendicular  to 
the  cross-flow  due  to  the  presence  of  the  jet  and 
afterbody  as  modeled  by  the  foil. 

i  *  1/4  the  distance  from  the  forward  focus  of  the 
foil  (corresponding  to  the  point  (-s.,0)  in  the 
C-plane)  to  the  trailing  edge. 

a  *  the  radius  of  the  transformed  circle.  This  term 
is  a  function  of  r,  the  distance  from  the  forward 
focus  of  the  foil  to  the  leading  edge. 


Appendix  A  contains  a  complete  derivation  of  the  above 


transformation . 

The  thickness  of  the  foil  is  a  function  of  both  r  and  i. 
As  described  above,  the  foil  length  is  a  function  of  the 
velocity  ratio  so  that  the  length  l  is  given  by: 


Initially  it  was  assumed  that  the  forward  focus  of  the 
foil  should  be  at  the  center  of  the  jet.  This  makes  r  equal 
to  one-half  the  jet  diameter.  When  substituted  for  the 
doublet  in  the  fully  corrected  model,  this  method  provided 
significant  improvement  in  the  downstream  area.  The  up¬ 
stream  correlation  was  poor,  however,  since  the  foil  was  too 
thick  and  too  blunt  at  the  leading  edge.  It  was  necessary 
to  adjust  the  value  of  r  to  provide  a  thinner  and  more 
streamlined  body.  Next  a  foil  with  its  thickness  equal  to 

one  jet  diameter  was  tested.  Using  the  formula,  maximum 

2 

thickness/d  *  5.2(m/2,  -  .  7(m/i)  )  [13],  where  m  is  a  function 
of  r.  It  was  found  that  if  r/d  *  . 12d/£,  where  i  is  given  by 
equation  (10) ,  the  thickness  is  nearly  equal  to  the  jet 
diameter  for  the  entire  range  of  velocity  ratios.  The  foil 
is  then  moved  forward  so  that  its  leading  edge  coincides 
with  that  of  the  jet.  This  configuration  (App.  B)  produced 
the  results  shown  in  Figs.  14,  15,  and  16  in  comparison  with 
experimental  values  [5]  for  velocity  ratios  of  4.0,  8.0,  and 
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10.0  respectively.  Radial  and  circumferential  pressure 
distributions  for  velocity  ratios  of  3.9  and  8.0  are  plotted 
against  the  experimental  findings  of  Feam  and  Weston  [5] 
in  Figs.  C-l  through  052. 

Attempts  were  made  to  use  other  fictitious  shapes  such 
as  ellipses  and  ovals  as  well  as  distributed  source-sink 
combinations.  The  foil  shape,  however,  produced  the  best 
results.  So  the  simple  relationships  for  the  near-field 
blockage  effect  given  above  were  retained  in  the  final 
model . 

Reviewing  the  results,  it  is  seen  that  for  a  velocity 
ratio  of  8.0,  the  correlation  between  this  model  and  the 
experimental  data  is  very  close  except  in  the  area  of 
highly  turbulent  flow  immediately  aft  (within  two  diameters) 
of  the  jet.  This  is  a  region  of  uncertainty  for  all  velocity 
ratios  due  to  the  three-dimensional  flow  caused  by  entrain¬ 
ment  ,  vortex  shedding ,  and  other  highly-viscous  wake 
phenomena.  Since  this  area  is  small,  it  contributes  little 
to  the  overall  thruster/hull  interaction.  For  a  velocity 
ratio  of  3.9,  the  correlation  is  excellent  except  for  the 
mid-  and  far- field  regions  aft  of  the  jet  centerline. 

This  was  expected  since  it  was  determined  earlier  that 
the  blockage  effect  of  the  jet  has  little  effect  in  the 
far-field.  It  was  observed  that  if  r  is  such  that  it  causes 
the  foil  to  have  a  circular  shape  at  the  leading  edge  the 
correlation  in  the  area  forward  of  the  jet  is  improved. 
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However,  the  agreement  deteriorated  in  the  areas  to  the 
sides  and  aft  of  the  jet.  By  varying  r/d  as  a  function  of 
the  location  of  the  point  at  which  the  velocity  perturbations 
are  being  calculated,  the  shape  of  the  foil  can  be  altered 
so  that  it  is  nearly  circular  at  the  leading  edge  and  becomes 
thin  aft  of  the  jet  centerline.  This  configuration  worked 
better,  but  it  failed  to  significantly  improve  the  overall 
pressure  distribution  resulting  from  r/d  =  . 12d/£.  This 
observation  did  lead  to  the  idea  of  using  a  half-cylinder 
shape  for  the  jet  blockage.  That  is,  the  fluid  forward  of 
jet  "sees"  a  cylinder  but  once  past  the  jet  centerline 
the  influence  of  the  cylinder  is  rapidly  diminished.  The 
results  of  this  configuration  are  presented  in  Fig.  17  from 
which  it  can  be  seen  that  while  this  is  a  significant 
improvement  over  Waterman's  model,  it  is  not  as  accurate 
as  the  full  foil-shape  model.  It  is  concluded  that  further 
improvement  of  the  model  must  involve  a  more  sophisticated 
treatment  of  the  vortices  within  the  plume  as  well  as  the 
blockage  effect  of  the  jet. 
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IV.  CONCLUSIONS 


As  evidenced  by  Figs.  14,  15,  16,  and  C-l  through  C-52 
this  analysis  adequately  predicts  the  pressure  distribution 
on  a  flat  plate  from  which  a  jet  is  emitted  into  a  crossflow. 
Application  of  the  model  is  limited,  however,  to  velocity 
ratios  VR  *  3.9  or  greater.  This  is  due  to  the  method  by 
which  the  length  of  the  fictitious  afterbody  is  determined. 
Above  velocity  ratios  of  about  4,  this  length  decreases  for 
increasing  VR,  but  plots  of  the  coefficients  of  pressure 
along  the  crossflow  axis  similar  to  Fig.  12  indicate  that 
the  afterbody  length  should  decrease  also  for  values  of  VR 
less  than  about  4.  This  is  consistent  with  the  observations 
of  Feam  and  Weston.  They  reported  that  an  anomaly  in  the 
pressure  distribution  in  the  wake  occurs  at  a  velocity  ratio 
of  3.9.  The  implication  is  that  the  rate  of  pressure  recovery 
in  the  wake  region  is  a  minimum  for  a  velocity  ratio  of  about 
3.9  [51. 

No  attempt  has  been  made  to  accurately  predict  the  pres¬ 
sure  distribution  in  the  highly  turbulent  wake  of  the  jet. 

This  area,  having  a  width  of  about  one  jet  diameter,  extends 
downstream  from  the  jet  along  the  crossflow  axis.  The  pre¬ 
dicted  values  of  the  plate  pressures  reported  here  cannot  be 
relied  upon  since  most  of  this  area  is  contained  within  the 
fictitious  afterbody. 
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Figure  18  shows  the  results  of  this  analysis  plotted  as 
the  ratio  of  the  induced  side  force  to  the  jet  thrust/ 

-Y  /Y. ,  versus  the  velocity  ratio  VR.  These  results  are 
compared  with  NPL  pressure  measurements  on  a  flat  plate  and 
the  model  tests  of  English  [1]  and  Chislett  and  Bjorheden  [2]. 
The  slope  of  the  predicted  -Y  /Y.  curve  is  approximately  5.5, 

3  t 

whereas  that  of  the  experimental  data  is  about  5.3  for  both 
the  flat  plate  and  the  model  tests.  While  this  analysis 
predicts  somewhat  larger  values  for  the  overall  (integrated) 
interference  effects,  the  small  difference  between  these 
slopes  indicates  that  the  nature  of  the  relationship  between 
the  interference  effects  and  the  velocity  ratio  predicted  by 
the  analytical  model  is  correct.  Feam  and  Weston  did  not 
report  overall  interference  effects  so  the  correlation  between 
their  work  and  the  NPL  experiments  is  not  known.  Since  the 
assumptions  concerning  the  shape  and  size  of  the  afterbody 
were  based  on  Fearn  and  Weston's  data,  the  degree  to  which 
this  analysis  is  in  error  cannot  be  determined  from  Fig.  18 
alone.  Even  with  these  deficiencies,  it  is  felt  that  this 
analytical  model  adequately  predicts  the  pressure  distribu¬ 
tion  on  a  flat  plate  and  can  be  used  as  a  first  approximation 
for  design  purposes. 
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V.  RECOMMENDATIONS  FOR  FURTHER  STUDY 

In  Chapter  III  it  was  noted  that  afterbody  shapes  other 
than  a  symmetric  foil  have  been  tried  without  significant 
improvement.  This  does  not  mean  that  this  approach  showed 
no  promise  for  possible  improvement  in  the  model,  but  only 
that  the  actual  shapes  tested  were  not  the  best.  It  is  felt 
that  this  model  can  be  improved  by  changing  the  shape  used 
for  the  jet  and  afterbody. 

As  reported  earlier,  the  crossflow  fluid  particles  near 
the  plate  apparently  see  a  cylindrical  shape  when  approaching 
the  jet.  As  they  pass  the  jet,  some  particles  are  entrained, 
some  join  the  turbulent  flow  within  the  wake,  and,  further 
out,  they  pass  as  if  the  jet  and  wake  together  are  a  solid 
object.  Just  downstream  of  the  jet,  along  the  crossflow 
axis,  there  is  an  area  of  high  entrainment.  This  causes  a 
high  local  pressure  drop  out  to  about  one  jet  diameter  for 
VR  *  8.  Pressure  recovery  is  rapid  past  this  point.  This 
leads  to  the  idea  of  using  a  shape  such  as: 


This  shape  has  a  cylindrical  leading  edge,  the  maximum 
thickness  is  just  downstream  of  the  jet  center,  and  then  it 


narrows  quickly  to  a  ralativaly  thin  tail.  The  total  length 
is  the  sane  as  that  for  the  symmetric  foil. 

It  is  expected  that  this  shape  can  be  formed  using  the 
Joukowsky  transform  of  the  existing  model  by  varying  r  as  a 
function  of  the  angular  position  of  the  particular  point  at 
which  the  perturbation  velocities  are  being  calculated. 
Although  a  functional  relationship  which  will  work  over  a 
range  of  velocity  ratios  has  not  yet  been  found,  it  is  pos¬ 
sible  that  one  exists.  This  idea  might  also  be  pursued 
using  panel  (distributed  singularity]  methods. 

Another  area  of  uncertainty  during  the  course  of  this 

study  was  the  empirical  jet  trajectory  equation.  Reports  of 

experiments  to  determine  the  jet  axis  provide  a  wide  range 

of  possible  trajectories.  If  the  trajectory  equation  used 

in  this  model  is  accepted  as  correct,  the  trajectory  con- 
2 

stant  B  *  . 19VR  stxll  may  not  be  accurate  over  a  wide  range 
of  velocity  ratios.  In  particular,  it  is  suspected  that  it 
may  be  in  error  for  velocity  ratios  near  4  where  an  anomaly 
in  the  plate  pressure  distributions  occurs.  The  cause  of 
this  anomaly  is  not  clear,  but  it  is  suspected  that  this  is 
the  point  where  the  jet  momentum  is  so  low  compared  to  that 
of  the  crossflow  that  the  jet  reattaches  to  the  plate  within 
a  few  jet  diameters.  If  so,  the  jet  trajectory  would  be 
quite  different  than  that  which  is  currently  modeled. 

Finally,  the  entrainment  of  the  crossflow  fluid  into  the 
jet  is  implicity  accounted -for  in  the  formulation  of  the 
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horseshoe  vortices.  However  #  the  pressure  drop  associated 
with  the  jet  entraining  the  surrounding  fluid  is  not  explic¬ 
itly  determined.  This  model  predicts  a  zero  pressure  distri¬ 
bution  at  very  high  velocity  ratios  (zero  ship  speed) . 

Further  improvement  of  this  model  in  accounting  for  entrain¬ 
ment  is  needed  if  it  is  to  be  used  to  predict  the  performance 
of  bow-thrusters  in  ships  with  no  way  on. 
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Figure  1.  Hull  low  pressure  areas  at  various  speeds. 
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Figure  2.  Hull  suction  effects  vs.  R 
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Figure  6.  Comparison  of  experimental  findings.  [5] 
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Figure  9.  Comparison  of  pressure  coefficient  contours,  experiment 
vs.  fully  corrected  conforming  vortex  model,  VR=8. 


Figure  10.  Comparison  of  pressure  coefficient  contours,  experiment 
vs.  fully  corrected  conforming  vortex  model,  VR=4. 


Fully  corrected 


I 


VR  :  10.0 


Figure  12.  Plots  of  coefficients  of  pressure  along  the  free 
stream  axis  used  to  determine  the  length  of  the 
foil.  Data  from  [5]. 


Figure  14.  Comparison  of  pressure  coefficient  contours,  experiment 
vs.  present  analytical  model,  VR=4. 


Figure  16.  Comparison  of  pressure  coefficient  contours,  experiment 
vs.  present  analytical  model,  VR=10. 
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Figure  18.  Comparison  of  theoretical  values  of  -Y  /Y.  with 
experimental  findings.  [II  sc 
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APPENDIX  A 


DERIVATION  OF  THE  JOUKOWSKY  TRANSFORMATION 

2 

The  Joukowsky  transformation,  Z  =  c  +  &  /c,  is  used  to 
transform  a  foil  in  the  z-plane  to  a  cylinder  in  the  s-plane 
so  that  the  foil  can  be  treated  as  a  doublet  in  a  uniform 
flow.  The  terms  used  in  this  transformation  are: 

Z  -  the  complex  position  vector  of  a  point  in  the 
Z-plane  Z  *  x  +  iy 

5  -  The  complex  position  vector  of  a  point  in  the 
C-plane.  c  =  £  +  in 

i  -  A  positive  real  constant  equal  to  1/4  the  distance 
between  the  foci  of  the  foil, 
a  -  The  diameter  of  the  cylinder  in  the  c' -plane, 
m  -  The  distance  between  the  n  and  n '  axes.  This  is  a 
function  of  r. 

r  -  The  distance  from  the  forward  focus  to  the  leading 
edge  of  the  foil.  That  is  r  *  )c  -  2l\  where  -c  is 
the  abscissa  of  the  leading  edge. 
z' -plane  -  The  linearly  transformed  c-plane  used  to  give 
the  foil  thickness.  ?'  =  s'  +  in' 

From  the  Joukowsky  transformation  and  Fig.  13  it  is  seen 
that: 

c  *  .  5  ( Z  ±  yj  Z2  -  4i2  )  (Al) 
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m  *  .5  (X  -  L) 

C1  =  c  - (-m)  *  Z  +  m 


The  complex  potential  of  a  cylinder  in  a  uniform 
is  given  by: 

u>U')  »  U^U'  +  a  A') 


The  complex  velocity  is  then  found  as  follows: 


U«  +  u 


.  ,  _  du  _  dw  d?'  dc 

1V  az  “  arc '  ar  ay 


dw  _  dw  dc '  /  dZ 

az  ■  a?*  ar  /  ar 


then  substituting 


and  since  c'  -  a  =  ?  -  s. 


dm 

ai 


(A2) 

(A3) 

stream 

(A4) 

(A5) 


(A6) 
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APPENDIX  B 


COMPUTER  PROGRAM 


XXXXXXXXXXXX XXXXXXXXXXXXXXXX 
X  X 

X  THIS  PROGR  AM  SIMULATES  A  X 
X  FLUID  JET  INJECTED  I  NT  P  X 

x  a  uniform  crossflow  and  x 
x  further  determines  th=  x 

X  PRES  SURF  DISTRIBUTION  ON  X 
X  THE  SURFACE  rROM  WHICH  X 
X  THE  JET  IS  EJECT FP  X 

X  X 

XXXX XXXXXXXXXXXX XXXXXXXXXXXX 

CEFINITICN  OF  TERMS 

NR  AO  -  NUMBER  OF  RAOUL  POSITIONS  PM  SURFACE 

(INPUT) 

NANG  -  NUMBER  OF  ANGUtAQ  POSITIONS  PN  SURFACE 

(INPUT) 

NPOINT  -  NUMBER  OF  ARBITRARY  POINTS  ON  SURFACE 
<NRAO*NANG) 

UJ  -  JET  VELOCITY  (  TNPUT) 

UM  -  MAIN  STREAM  VELOCITY  (INPUT) 

O  -  ACTUAL  OM»FT«  OF  JET  npjcjfc  ( INPUT) 

R  -  NON-O  TMCNSI  PNA)  P  AOIAL  POSItTPN  OF  SUR  F  ACe 

POINTS  (DIMENSIONED  at  LEAST  NR  AO)  (INPUT) 

THETA  -  ANGULAR  POSITION  oc  SUP  FAC*  POINTS 
(DIMENSIONED  AT  LEAST  NAnq) 

XP  -  X-°0$  T  T  I  ON  OF  SURFACE  ®OTNTS 

(DIMENSIONED  AT  LEAST  NPOINT) 

YP  -  Y~°OS  ITI  ON  OF  SURFACE  POINTS 

(01  MENS  10 NEO  AT  LEAST  NPOINT) 

X  -  HORIZONTAL  COORDINATE  OF  JET  AXIS 

Z  -  VERTICAL  COORDINATE  OF  JET  AXIS 

ZINC  -  NON-OI^ENSIONAL  INCREMENTAL  STS®  SIZE  FOR  Z 
OELTAZ  -  Z-POSITION  OF  EFFECTI VE/VIRTUAL  SOURCc 
HKSAV  -  CUMULATIVE  VORTICITY  ALONG  THE  JeT 

XI  -  NATURAL  COORDINATE  TANGENT  TO  JET 

ZETA  -  NATURAL  COORDINATE  NORMAL  TO  JeT 

ETA  -  NATURAL  COORDINATE  PERPEND! CULAR  T0  XI 
AND  ZETA 


AL  °H  A 


ANGLE  BETWEEN  X-AXIS  AND  XI-AXTS 


nnononDnnnoonnnnn 


PA R 5 A  -  INCREMENTAL  PLATS  AREA 

TOTFCR  -  TOTAL  RESULTANT  FORCE  ON  PLAT- 

XRAR  -  X-POSITION  OF  RESULTANT  FQ3Cp 

YT  -NORMALIZED  INDUCED  FORCE 

ELL  -  1/4  THE  DISTANCE  BETWEEN  THE  ROC! 

OF  THE  FOIL  USED  FOR  THE  BLOCKAGE  EFFECT 

RA04  -  THE  DISTANCE  FROM  THE  FORWARD  FOCUS  OF  THF 
FOIL  TO  THF  LEADING  EDGE 

ZCP  -  THE  COMPLEX  VELOCITY  INDUCED  IN  The  UNIFORM 
FLOW  OUE  TO  THE  RRESENCS  OF  THF  FOIL 


IMPLICIT  REALMS  I A-H 

oi  mansion  uc  so  .so)  ,v 

X  VFL (450) ,CP(50,50) 
OIMFNSICM  XP (450 ) * Yd 
COMPLEX  ZED* 16 • ZST*L 
*ZE  l*  16  »Z4*  16  *  ZH4*L 6, 
DATA  INPUT 

RE  AD ( 5  *60 )  NPAO 
R E A 0  ( 5  »  60  )  NANG 
6 C  FORM  AT  ( T4) 

NP 01  NT  *  NRAD*NANG 
ZERO  OUT  STORAGF  MATRI 
DO  30  I  *  1 ,  NP.  AD 
'00  20  J  *  1 ,  NANG 
U(  I  *  J)  *  0.0 
V(  I.  J)  *  0.0 
c° (I  .J  )  *  n.0 

C3NT  TMUe 
CONTINUE 

RE AO (5. 40  I  UJ.UM.D 
WRITE(6,50J  UJ,UM,0 
VELOCITY  RATTO 
VR  * (UJ/UM) 

RE  AD (5. 70}  ZINC 
RE  AO ( 5  »  70 ) ( R ( I ) ,1*1, 
FORMAT  (F10.5  ) 

SET  CONSTANTS  UScO  IN  F 
Bl  *1.0 
BHA  *  0.5 
84  *  4.C 

SET  POIL  "ARAMETERS 
ELL  *  6.0/VR 
RAO  J  *  .12  /ELL 
ALAM  »  .5*(2.*ELL*RA 
*4. 0* ELL** 2)1 
EM  *  0.5* (ALAM  -  ELL 
POSITION  OF  EFFECTIVE 
0F  L  T  A  Z  *  .6*0SQRT(VR 
WRITF(6,a0»  DFLTAZ 
WRITE(6,90J  ZINC 
FO  RM  AT  (  3F  3 .4  ) 

FQRMAT( ’O1 ,2X, *  JET  v 
X»VELCCITY  *  »,F3.4,2 


,0-Z) 

(50,50  ) ,  W(  50  ,50 )  ,R(  50),  THETA  (50  ) 

(450) ,  Z° (450  ) 

6, ZC°* 16 ,ZSTP*l6 , ZR*1 6, 

ZM*I6 


NRAD) 

CIL  BLOCKAGE  CALCULATIONS 


FORMAT! 


!Ml£ 


X’DfAMFTERS  ABOV*  THC 
FORMAT  ('0»  ,2X  , 1  Z  ST 
DFFINITIDN  of 

PI  *  2.0*ARS IN ( 1 ,000 
trajectory  constant 

3  *  0.19*VR**2 


DJ+OSO  RT ( ( 2.*ELL+R AOJ )**2- 


SOURCE 

) 


ELDCTTY  *  »,  FB. 4,2X, 'UNIPO^Mt , 
X,»JET  DIAMETER  *  *,F3.4) 

TIVE  SOURCE  I  S' , F5 . 2 , 2X , 

PL  ■*  TE  • ) 

EP  SIZE  *  *  ,  c5 . 3  ) 

00000) 


CALCULATION  OF  INTERFE0  FNCE  VELOCITY 
L  *  0 

CALCON  x  OFLOAT(  NANG  )  -  t.O 
INCREMENT  RADIUS  (INPUT  VALUES) 

DO  140  I  *  l,  NR  AO 
PC  AL  C  =  ©(I) 

INCREMENT  THETA  (5  DEGREE  INCREMENTS) 
00  130  J  »  1  ,  NANG 

K  x  j-l 

THETA(J)  «  DFLOAT(K)*  (oj/CALCCN) 
ARAD  x  ALAM  -  EM 
SHIFT  TO  CARTESIAN  COORDINATES 
L  x  L+l 

XP  (  L  )  *  RCALC*OCOS  (THET  A(  J  )  ) 

Y°  (  L  )  *  RCALC*OSIN<THETA(  J  )  ) 

HKSAV  *  0.0 
IZ  =  0.00001 
100  21  »  Z2 

Z2  *  Z1  +  ZINC 
Zll  «  Z1  +  DELTAZ 
Z22  »  Z2  ♦  OELTAZ 

XI  a  B*(DCnSH(  Zl/8)-1.3) 

X2  x  B*(OCQSH( Z2/8I- 1.0 ) 

ALFONF  a  OATANI 1.0/0  StNH( Z l/B)  ) 

AL FT WO  *  0ATAN(1.0/0SINH(Z2/9)  » 
ALFBAR  x  OATAN(  (  Z2-Z  l)/ (X2-X1)  ) 

SINAI  x  OSIN(ALFONF) 

SINA2  x  OSIN(ALpTWO) 

SINAB  x  OS  IN  (  ALC  BA  R ) 

C0SA1  *  nCOS ( A LcONE) 

C0SA2  =  OCOS(ALcTWO) 

CO  SAB  *  OCOS(ALERAR) 

NATURAL  COOP  0  TNATFS  IN  TERMS  OE  CART=S' 
XI  ONE  x  -  <  Zli*S!NAl  ♦  (  XI-  XP  ( L  >  )  *C0! 
XI  TWO  *  -  (  Z22*  S  1  N A2  ♦  (  X?- X*»  (1  )  )  *r.r>j 
XI8API  *  ~(Z1.1*SINAR  (U-XP(L)HT/ 
XTRAR2  =  — ( Z22  IN AB  +  ( X2 -XP ( L )  ) *C' 


TERMS  OF  CART=SIAN  COORDINATES 
+  (  Xl-XPU.  )  JvCOSAl  ) 

+  (  XZ-X'Ml  )  )*r.ns  \Z  ) 

(n-xpumcosAB) 

+  (X2-XP( L ) )*CDSAB) 


ETA  x  YP(  L  ) 


ZETA1  x  (Zll*COSAl 
ZETA2  =  (Z22*COSA2 
ZETA9  »  (Zll*C0$4B 
GEOMETRIC  CONSTANTS 
ET  PL  *  =TA  ¥  0.5 
ET Ml  a  eta  -  0.5 
SA  x  X  TONE** 2 
SB  *  XITW0**2 
SC  x  xi£*ARl**2 

50  *  X  I  BAR  ?**2 
SE  *  Z^TAl**? 

SF  x  ZCTA2**2 
SG  *  ZETAB**2 
SH  *  ETPt.*+2 

51  a  PTMT**2 

OA  a  OSQPTISA  ¥  Sc 
OB  *  OSCR T{ S A  +  SE 
OC  x  OSQP  T(  SB  *  SE 
OE  *  DSQR  T ( S  B  ♦  SE 
OF  *  0  S OR  T  (  S 0  ♦  SG 


(Xl-X  P(L)  ) *S  TNA I ) 
(  X2-X  P(  L)  )*SINA2) 
JXl-XP(L)  )*SIN4B) 
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DC 

9 

DSCRKSC  ♦ 

SG 

DH 

9 

OSCPKSO  ♦ 

SG 

01 

9 

OSCRT (SC  ♦ 

SG 

CONI 

■  (ETOL/OA 

— 

CON  2 

«  (FTPL/OC 

— 

COM3 

a  ( X I  BAR  2/DF 

C0N4 

»  (XI3AP2/0H 

Ul 

9 

ZETA  i*CON  l 

U2 

9 

ZFTA2*CON2 

V3 

9 

ZPTA8*C0N3 

V4 

9 

ZETA8*C0N4 

wi 

9 

X  ICNF*CON 1 

U2 

M 

XI TuO*C ON 2 

W3 

M 

ETMI*r.0N3 

W4 

3 

ETRL*C0N4 

SI  ) 
SH  ) 
SH) 


FT  M  I/OR  )  /(  SA  «■  SE> 
ETMI/OE)  /  (SB  ♦  SP) 

-  xib  ari  nr,)  /(  sg  sn 

-  XI0 A3  i  /on/  (  SG  ♦  SH) 


INCREMENTAL  VORTICITY  TERM 

FF  *  (  |0EXP(Z2/B)-DFXP(  Zl/B)  )/{  l.OOSXPt  (Z2+ZU/8)  )  ) 
HK  *  (1. 0/8.0)*  VR**2  *  QATAN(PF) 

HK  »  HK  ♦  HKSAV 
HK  SA V  *  HK 

INCREMENTAL  INTERFERENCE  VELOCITY  COMPONENTS 
UP  *  0.0 

UP  *  HK*<  (-WI*SINAH-W2*SINA2*(-w3+k4)*SIMAB+UI*CPSU) 
X  - ( U2*CPSA2 ) ) 

VP  *  0.0 

VP  3  HK*(V3  -  V4 ) 

U(  I  *  J )  =  U ( I  VJ  I  ♦  UP  *2.0 

v(  r,jj  *  vti.j )  vp *2.0 

TEST  FTP  CON VERGFNCE 

TEST  =  DABS(UP)  ♦  PABS(VR) 

IF  (TEST. LE.O. 000001 )  GQ  TO 
GD  TO  ICO 
CO  NT  IN  Uc 


110 


110 

C  CALCULATION  OF  8LOCKAGF  VELOC  ITY  TERM 
C  TRANSFORM  COORDINA^S  TO  Z  »L  ANE 
83  AZ  =  XP(L)-  2.0<:ELl  ♦  (.S-RAOJ) 

IF  (AZ.EC.o.O.ANO.YR(  L).cQ.0.0)GD  TO  1L5 
C  CHANGE  VARIABLES  TO  COMPLEX  NUMBERS 

az  a  yp(l> 

aa  a  o.o 

ZED  a  QC.MPLX  ( A Z  *  8Z  ) 

ZR  a  OCMPL  X  (  ARAO  *  flB) 

ZFL  *  OCMPI.X  (ELI  *  8B) 

Z4  a  0C«PLX(B4,9B) 

Z«  «  OC  yPLX(  EM  *  88  J 
ZHA  »  OCMPLX  (BHA »  88) 

IF(4Z.Lt.0.0  )G0  TO  89 

C  TRANSFORM  POSITION  VECTOR  IN  Z  PLANE  TO  Z  FT  A  PLANE 
ZST  «  ZHA*(ZEO  +  COS QRT ( ZE D**?-Z4*ZEL**2 ) ) 

GO  TC  91 

89  ZST  «  ZHA*( ZED  -  CDS QRT( ZED**2-Z4*ELL**2) ) 

C  TRAMS  FORM  POSITION  VECTORS  FORM  ZSTA  P L A NF  TO 
C  ZETA  PRIME  PLANE 


91  ZS TR  a  ZST  ♦  ZM 
SET  VELOCITIES  TO  ZERO 
XEO  «  ORE  AL  (  ZST°  ) 

YEP  ■  PIMAG(ZSTP) 

PR  a  0S0RT(XEP**2  ♦ 
IF(PR.LF.(  ARAO/14.0) 


WHEN  WITHIN  THE  FOIL 


YE0**2  ) 
JGO.TO 


95 


ir  t  L  1  ItlUI  I  _ 

CALCULATE  the  DERIVATIVE  OF  THE  COMPLEX  POTENTIAL 
95  ZCP«((Z$T/ZSTP)**2)*(ZST°  ♦  ARAOJ/IZST  +  RLL ) 

EXTRACT  THF  HOWIZ.  ANO  VFR  T .  PE  RT  Aft AT  1 0  N  VCL^CITIES 
UC  »  0REAU7CP)  -  1.0 
VC  a  -OIMAG(ZCP) 

GO  TC  il‘> 
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115  tC  *  O.C 
VC  »  0.0 

119  u<  r. j)  *  u(  i  *j )  +  uc 

VI  I  .J)  *  VI  T  ,J)  ♦  VC 
C  CALC'JLATICN  OF  PRESSURE  CQEFF  IC IENT 

120  C  °(  I  •  J  )  -  -2.0*(U(I  ,JI  +  ,5*il  C  t  *  J  1**2  ♦  .5*V(  I ,  J  )**2) 
C  CHANGF  DESIGNATION  OF  VARIABLE 

ZP  (l  )  *  CP(  I.J  ) 

VEUL)  *  OSORTIIUI  I«  JJ+1.0  )**2  +  V(I,J)**2) 

130  CONTINUE 
140  continuf 

C  CALCULATION  OF  INCREMENTAL  FORCE 
TOTFOR  *  0.0 
TYMOM  *  0.0 
NAREA  x  MRAD  -  1 
TO  170  1*1,  NAREA 
K  x  I  ♦  l 

pi  =  r  <  n 

PO  x  RIK) 

C  CALC'JLATtCN  OF  INCREMENTAL  AREA 
AREA  *  PI*(PQ**2  -  R  1**2 ) 

P  AREA  =  AP?A/(2.0*CALC0N) 

DO  160  J* l, NANG 

C  C  ALCULA  TJ  CM  OF  AVERAGE  PRESSURE  COEFFICIENT 
CPAVG  *  (CPU.J)  +  CPIKt J  )  )/2.0 
C  CALCULATION  OF  FF  FEC.T  IV  E  FORCE 
FORCE  *  P  APE  A* CPAVG 
IF (J.FO.l I  GO  TO  150 
IF (J. EC. NANG)  GO  TO  150 
FORCE  =  2  •  0*FQRCF 
C  CALCUIATICN  OF  MOMENT  ARM 
150  XM  =  (IPG  +  RI )/2.0) *0COS( TMFTA(J) ) 

C  CALCULATION  OF  MOMENT  ABOUT  Y-AXIS 
vmom  a  FORCF#XM 
C  SUMMATION  OF  TOTAL  PQRC  E 

TOTFOR  *  TOTFOR  ♦  FORCE 
C  SUMMATION  CF  total  MOMENT 
TYMOM  =  TYMOM  +  YMOM 
160  CONTINUF 
170  CONTINUE 

C  CALCULATION  OE  CENTER  CF  ACTION 
X3AR  x  TYMCM/TQTFOR 

C  CALCULATION  OF  NORMALIZED  INDUCED  FORCE:  Y/T 
VRINV  x  1.0/  VR 
YT  x  2.0*TOTFOR*VRINV**2/P  I 
C  OATA  OUTPUT 

WR I T  EC  6.230 )  TOTFOR*  X5AR 
wR  ITE(  6  *240)  YT,  VRINV 
IdO  00  2  20  I  *  l « NANG 

ANGIE  x  THCT A(  I  )  *1  80  .O/PI 
WR I T E I  6  *  1 9 0 )  ANGLE 

190  c0  PMAT  ( 'O' ,2X» 'THETA  *  • ,F 6. 2. 8X, • R * , 9 X, • CP* *9Y , ' U* , 

X  9X,  'V»  ) 

DO  210  J  *  1 ,NPAD 

WR ITEI 6*200)  R{J),CP(J,I) ,UI J, I) *V(J,I) 

200  FORMAT!'  •  ,  23X  ,  F6 .2 ,  2X,  F8 . 2, 3X  ,  F  6 . 2,  4X  ,  F  6 . 2) 

210  CONTINUE 
220  CONTINUE 

230  FORMAT! 'O' ,2X, 'RESULTANT  FORCE  OF  ' , IX, *7.2, IX, 

X'ACTS  AT  X  =  '  ,F6.2) 

240  FORMAT { 'O' » 2X, 'Y/T  =  » , F6 . 2 , IX , • FOR  UM/UJ  »  ',F6.3I 
WRITE!  7*2  50)  (XPl  I  )  ,YP(I  ),Z°(  I)  ,l=i,NPOINT) 

250  FORMAT! 2F20. 4) 

270  STOP 
END 
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COMPARISON  OP  PRESENT  MODEL  WITH 
EXPERIMENTAL  RESULTS 


RADIAL  DISTANCE  CX/D) 


STANCE  (X/D) 


aynssayd  do  1N3I3I3J3Q3 


radial  lines,  VR  =  8,  theta  =  20  degrees 


Figure  C-4.  C  's  along  radial  lines,  VR  =  3,  theta  =  30  degrees. 


STANCE  (X/D) 


Figure  C--o.  C  *3  along  radial  lines,  VR  =  8,  theta  =  50  degrees 
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Figure  C-8.  C  's  along  radial  lines,  VR  =  8,  theta  =  90  degrees 
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Pigure  C-10.  C  's  along  radial  lines,  VR  *  8,  theta  =  120  degrees 


Figure  C-ll.  C  's  along  radial  lines,  VR  =  8,  theta  =  135  degrees 
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RADIAL  DISTANCE  (X/D) 


Figure  C-14.  C  ’s  along  radial  lines,  VR  =  8,  theta  =  170  degrees 


radial  lines,  VR  =  8,  theta  =  180  degrees 


Figure  C-16.  C  ' s  along  circumfrential  arcs,  VR  =  8,  radius  =  .58  dia 


Figure  c-17.  C  's  along  circumfrential  arcs,  VR  =  8,  radius  =  .75  dia. 


Figure  C-19.  C  ’s  along  circurafrential  arcs,  VR  =  8,  radius  *  1.4  dia 


Figure  C-20.  C  1 s  along  circumf rential  arcs,  VR  =  8,  radius  =  1.6  dia. 


Figure  C-21.  C  ' s  along  circuraf rential  arcs,  VR  =  8,  radius  =  2  dia. 


Figure  C-22.  C  ' s  along  circumf rential  arcs,  VR  »  8,  radius 


circumfrential  arcs,  VR  =  8,  radius  =  3  dia 


circumfrential  arcs,  VR  =  8,  radius  3.5  dia 


Figure  C-25.  C  ' s  along  circumfrential  arcs,  VR  =  8,  radius  =  4  dia 


Figure  C-26.  C  ' s  along  circurafrential  arcs,  VR  =  8,  radius 


radial  lines,  VR  =  3.9,  theta  =  0  degrees 
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Figure  C-30.  C  ' s  along  radial  lines,  VR  =  3.9,  theta  =  30  degrees 
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Pigure  C-31.  C  's  along  radial  lines,  VR  =  3.9,  theta  =  45  degrees. 


STANCE  (X/D) 


Figure  C-33.  C  'a  along  radial  lines,  VR  =  3.9,  theta  =  75  degrees 
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radial  lines,  VR  =  3.9,  theta  =  135  degrees 


radial  lines,  VR  =  3.9,  theta  =  150  degrees 


radial  lines,  VR  *  3.9,  theta  ■  170  degrees 


Figure  041.  C  ' s  along  radial  lines,  VR  =  3.9,  theta  =  180  degrees 
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Figure  C-42.  C  's  along  crlcumfrentlal  arcs,  VR  =  3.9,  radius  »  .58  dla 


Figure  C-43.  C  's  along  circumfrential  arcs,  VR  =  3.9,  radius  **  .75  dia. 
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circumfrential  arcs,  VR  *  3.9,  radius  =  1  dia 


Figure  C-45.  C  's  along  circurafrential  arcs,  VR  =  3.9,  radius  =  1.4  dia 


circumf rential  arcs,  VR  =  3.9,  radius  =  1.6  dia 


Figure  C-47.  C  's  along  circumfrential  arcs,  VR  =  3.9,  radius  =  2  dia. 


ircumfrential  arcs,  VR  =  3.9,  radius 


Pigure  C-49.  C  's  along  circurafrential  arcs,  VR  *  3.9,  radius  =  3  dia 


Figure  C-50  C  *s  along  circumfrential  arcs,  VR  =  3.9,  radius 


Figure  C-51.  C  'a  along  circurafrential  arcs,  VR  =  3.9,  radius  =  4  dia 


C-52.  C  'a  along  circurafrential  arcs,  VR  =  3.9,  radius 
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